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There has been substantial success in predicting the heats of formation of binary alloys both with ab-initio 
and semiempirical methods (1-8). There is still a good deal of progress to be made, however, in that 
experimental trends in the cohesive energies, heats of formation and lattice parameters as a function of 
composition are often not accurately reproduced. 

Recently, a new method for the computation of alloy formation and defect energies was introduced by 
Bozzolo, Ferrante and Smith (BFS) (7). We tested the method ability to predict the heat of formation as a 
function of alloy composition for disordered binary alloys. In this paper we present a further property of 
interest, the variation of the lattice parameter with composition. At present, Vegard's law (9), which is a 
simple compositional linear average, is the primary method of choice. However, it is well known that the lattice 
parameter variation often deviates from linearity (10). In this paper we present the predictions of the BFS 
method for the lattice parameter. A second goal of this paper is to provide a simple formulation of the 
algorithm for performing defect calculations. We will see that accurate prediction of the composition dependence 
of lattice parameters can be made with our simple method. 

BFS (7) requires information only about pure metal properties and certain experimentally determined alloys 
properties. We build on the formulation of equivalent crystal theory (ECT) (11-12), by dividing the total 
binding energy of the alloy into a chemical energy and a strain or structural energy (see, e.g., Refs. 2-4). We 
now outline the procedure for calculation of the heat of formation and lattice parameter versus concentration 
for fee binary alloys. We first apply the method to ordered alloy structures and then use the Connolly-Williams 
(CWM) (13) procedure for relating the ordered compounds to the disordered ones. 

Consider a unit cell containing N x atoms of atomic species X, (X = A,B,...), so that the total number of 
atoms in this cell is given by N = E x N x . The heat of formation of this cell is 

A^cell = E cell “ ]C N x E x (j) 

where E ce „ is the total energy of the cell and E x is the energy of an atom of species X in a pure crystal of 
its own species. IF E'(i,X) denotes the energy of the i th atom in the cell (i = 1,...,N X ) of species X then 

^i, = EEh‘W-Ex = EEu ( 2 ) 

X i~l X i=l 

In our simulation of the process of alloy formation, each individual atom undergoes two separate 
transformations with respect to its initial state given by a perfect, pure crystal of its own species. Each 
transformation contributes to the heat of formation of the alloy. 

The first process (strain energy) deals with the structural change arising from neighbor locations in the alloy 
being different from in the elemental single-crystal environment. In order to compute the strain energy of the 

l-th atom of species X, e. ^ we just ‘flip’ every atom in the actual alloy structure, surrounding atom (i,X), into 
an atom of the same species X, and perform a regular ECT calculation (see Ref. 12 for details). Then, 
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where a e denotes the equilibrium lattice constant of a pure X crystal, E c is the cohesive energy and l x is a 
scaling length (12). The equivalent lattice parameter a^ x is determined by solving the appropriate ECT 
equation applied to atom (i,X), in the defect (but pure) crystal (12). For the simple case considered in this 
work, a. x is just the actual alloy lattice parameter. The second process deals with the change in chemical 
composition. The chemical energy contribution to the heat of formation from atom (i,X), eP x , is the energy 
change due to some neighbors of atom (i,X) being of a different atomic species. Here, we ‘freeze’ the neighbor 
locations to be those corresponding to a pure X crystal allowing only for a change of identity of the neighboring 
atoms. Rigourously, this change in chemical composition cannot be considered a ‘defect’ in the context of ECT. 
However, to a good approximation, the underlying concepts of ECT should remain valid in the case of alloys, 
and we will adopt a similar formulation to describe the perturbation due to the dissimilar atomic species. We 
assume that the global property parameterized by the ECT parameter or (12) (i.e., the tails of the overlapping 
electron densities) can be separated into pairs of interacting atoms. In this approximation, the electron density 
in the region between two atoms of the same species would be unaffected by the presence of neighboring atoms 
of a different species so that the perturbation would be localized in the region between two dissimilar atoms. 
This assumption enables us to define the parameter <*xy as a XY “ a X + ^YX> w h ere c*x a va ^ ue 

for the pure metal X and Ay^ is a correction introduced by the presence of a neighbor of species Y. The 
parameters Ayx an ^ A X y are ^ e on ^Y new parameters introduced in our method, and they will be 
determined by fitting to appropriate experimental data. 


With the introduction of the parameters Ay X and A X y we can compute the chemical energy of atom (i,X) 


using ECT. The chemical energy e i,x is obtained by performing two similar ECT calculations: 


c 

e. v — e 
i,X 


£x( A Y,x) 


< : x<0). 


( 4 ) 


For the first term, as explained above, the chemical composition is included in the appropriate value of Ay x 

C / A \ „X , /, CM -a^x c* / C X\ ( 5 ) 

e i,x( A Y,x) _ li.xEcl 1 r + a >.xj e f a i,X " \ a i,x a e J / X 

c* • * c . . 

with 7 j x = 1 if a j x > 0 and 7j x = 1 otherwise. The equivalent lattice parameter a. x is obtained by solving 

the following ECT equation (12) 


-|o x + — 1 

„ n P x -a v R. I A v I 

NR ( e x 1 + M H 2 e ' x ' 

-U YX ♦ JL r2 

Nvvr Px e “ Yxr ' + Mwr!* e ^ ^ 


( 6 ) 


£ l, YX r i 
Y 


where N (M) is the number of nearest-(next-nearest)-neighbors in a perfect crystal of species X,R 1 (R 2 ) is the 
nearest-(next-nearest)-neighbor distance in the equivalent crystal of lattice parameter a C y . The sum on the 
r.h.s. of Eq. (6) is over the atomic species Y of the Ny X (Myx) nearest- (next-nearest) -neighbors of atom (i,X) 
located at a distance r^rj). The second term in Eq. (4), ej^O) involves a calculation similar to Eqs. (5-6) but 
with Ay X = 0. This is done in order to free the chemical energy from any structural defect (e.g., a surface) 
thus retaining only the contribution of the chemical composition of the surroundings of atom (i,X). 
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Finally, the contribution 


e i,x 


of atom (i,X) to the heat of formation (Eq. (2)) is 


e i,X 


*,X 


gi.X 4 


i.X 


( 7 ) 


where the coupling factor between the strain and chemical energy contributions, gj x “ e a ' ,x ensures that the 
chemical energy vanishes at large interatomic distances. 


In this work, we are concerned with the fcc-based disordered binary alloys A x B 1 _ x , which, in keeping with 
CWM, we will compute from the corresponding ordered structures A m B 4 m (m = 0(B,fcc); l(AB 3 ,Ll 2 ); 
2(AB,L1 0 ); 3(A 3 B,Ll 2 ); 4(A,fcc)). If no relaxation of the individual atomic sites is allowed, then the strain 
energy, e x , is given by Eq. 3 with a x = r, where r is the actual interatomic distance. Within this 
approximation, the second term in the chemical energy (Eq. (4)) vanishes, since there are no structural defects. 
For a given ordered structure m, the ECT equation for the equivalent lattice parameter a C ^ (Eq. (5)) is 


obtained from Eq.(6) with M=12, M— 6, Rj = ^ R 2 ; R 


v r i ~ 


a 


r 0 and r 0 


The parameters p^, 


« x , and A x are listed in Ref. 12 and the coefficients N xx , N XY , M xx , and M XY are obtained from the 
number of nearest and next-nearest neighbors of each species for all the possible structures A m B 4 . The excess 
energy, AE m (r), as a function of the lattice parameter r is easily computed with Eq. (2) and, following 

Connolly- Williams approach (13), the heat of formation is obtained from 

AE D (r,x) = £ Wx m (l - x)<- m AE m (r) (8) 

m V J 

by finding, for each concentration x, the value of r that minimizes AE D (r,x). For different alloys A - B, the 
parameters A ab and A 0a were determined by reproducing the experimental heats of solution E BA and 
E AB- In Tnble I, we list the parameters A- for some binary alloys of Ag, Ni, Cu, Au, Pd, Pt, and Al, 


Table I. Experimental values of the heats of solution E AB and E 0A used to determine 
the parameters A ab and A ba for some binary alloys of Al, Cu, Ni, Ag, Au, Pd, and Pt. 


A 

B 

e ab 

e ba 

a ab 

a ba 

A 

B 

e ab 

e ba 

^AB 

^BA 

Ni 

Pd 

-0.088 

0.057 

-0.0401 

-0.04665 

Ag 

Pd 

-0.108 

-0.289 

-0.0431 

-0.02033 

Cu 

Ni 

0.09 

0.03 

-0.0131 

0.02395 

Cu 

Pt 

-0.299 

-0.532 

-0.0568 

-0.0444 

Cu 

Pd 

-0.392 

-0.436 

-0.04205 

-0.04795 

Ag 

Au 

-0.16 

-0.19 

-0.0313 

-0.0219 

Au 

Pd 

-0.195 

-0.355 

-0.0439 

-0.0348 

Cu 

Ag 

0.25 

0.39 

-0.0321 

-0.0394 

Ni 

Pt 

-0.330 

-0.282 

-0.0603 

-0.0529 

Cu 

Au 

-0.126 

-0.19 

-0.0588 

-0.05095 

Cu 

Al 

-0.20 

-0.35 

-0.0626 

-0.0526 

Ni 

Au 

0.22 

0.28 

-0.0614 

-0.0512 


as well as the experimental heats of solution used in this work. The experimental heats of solution are estimated 
from the experimental heats of formation curves in Ref. 14. A complete list of the parameters needed to 
calculate defect energies once the geometry and composition are known, and can be found in Ref. 12. In order 
to apply the method for a specific defect, e.g., a surface, one needs only know the position of each atom and the 
nature of their neighbors (i.e., A or B). 
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In Fig. 1 we show results for Ni-Cu, Cu-Pd, Ni-Pd chosen because of the different behavior exhibited by 
each one of these compounds. We display the heat of formation, compared to experiment (14), and the lattice 
parameter obtained with our approach. Figure 2 shows our predictions for the lattice parameter of other binary 
systems. We also compare these results with experiment (10) and Vegard's law (9). Note that our method 
is able to reproduce the particular features of the Ag-Au system, where, for a certain range of concentrations, 
the lattice parameter shows a contraction from the equilibrium lattice parameter of the pure components (10). 
We reproduce this exceptional behavior rather accurately, which is surprising since we are dealing with 
differences of hundredths of angstroms. In all cases, the agreement with experiment is excellent. In Figs. 1(a) 
and 1(b) we also present a comparison with an appropriate version of the embedded atom method (EAM) (4,16), 
Miedema's empirical approach (15) and experiment (14), for Ni-Pd. There is another EAM calculation of Ni- 
Pd (5). However, the parameterization used there is different from the original EAM formulation for alloys (8). 

In conclusion, we have developed a new semiempirical procedure for the concentration dependence of the 
heats of formation and lattice parameters of binary alloys consistent with ECT. This method accurately predicts 
the experimental behavior qualitatively and quantitatively. 
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Figure 1 .—{Left column) Comparison between the heat of formation as 
a function of concentration for Ni-Pd f Cu-Ni and Cu-Pd alloys: the 
solid curve indicates the results obtained in this work and the solid 
squares indicate the experimental values. The Ni-Pd figure also 
displays the results obtained with Miedema’s approach {dashed 
curve) and Johnson’s EAM (4,1 6) results (dash-dot curve). (Right 
column) comparison of lattice parameters of Ni-Pd, Cu-Ni and Cu-Pd 
alloys as obtained in this work (solid line) and the corresponding 
experimental values (solid squares). The Ni-Pd figure also displays 
the results obtained with Miedema's approach (dashed curve) and 
Johnson's EAM (4,16) results (dash-dot curve). 
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Figure 2. — Comparison between experimental values (solid squares) and the 
results obtained with this approach (solid lines) for the lattice parameter of 
several binary alloys of Ag, Ni Al, Au, Pd, Cu, and Pt. The results predicted 
by Vegard's law (dashed line) are also shown. 
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